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Abstract. Characterizing the shape and evolution of pul- 
sar radio emission beams is important for understanding 
the observed emission. The various attempts by earlier 
workers investigating beam shapes have resulted in widely 
different conclusions. Using a carefully selected subset of 
the recently published multifrequency polarimetry obser- 
vations of 300 radio pulsars (Gould & Lyne, 1998), we 
attempt to model the shape of pulsar beams. Assuming 
that the beam shape is elliptical, in general, and that it 
may depend on the angle between the rotation and the 
magnetic axes, we seek a consistent model where we also 
solve for the dependence of the beam size on frequency. 
From the six-frequency data on conal triple and multiple 
component profiles, we show that a) the pulsar emission 
beams follow a nested cone structure with at least three 
distinct cones, although only one or more of the cones may 
be active in a given pulsar; b) each emission cone is illu- 
minated in the form of an annular ring of width typically 
about 20% of the cone radius. 

Although some slight preference is evident for a 
model where the beam is circular for an aligned rotator 
& latitudinally compressed for an orthogonal rotator, 
the possibility that the beam shape is circular at all 
inclinations is found to be equally consistent with the 
data. While the overall size scales as p-O-^ (where P is 
the pulsar period) as expected from the notion of dipolar 
open field lines, we see no evidence in support of the 
beam shape evolution with pulsar period. 



Key words: pulsars: general-emission mechanism 



1. Introduction 

Most widely accepted emission models assume that pulsar 
radiation is emitted over a (hollow) cone centered around 
the magnetic dipole axis. The observed emission is gener- 
ally highly linearly polarized with a systematic rotation of 
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the position angle across the pulse profile. This behaviour, 
following Radhakrishnan & Cooke (1969), is interpreted in 
terms of the radiation being along the cone of the dipo- 
lar open field-lines emerging from the polar cap, and the 
plane of the linear polarization is that containing the field 
line associated with the emission received at a given in- 
stant. During each rotation of the star, the emission beam 
crosses the observers line-of-sight resulting in a pulse of 
emission. The observed pulse profile thus corresponds to a 
thin cut across the beam at a fixed rotational latitude. The 
information on the beam shape as a function of latitude, 
although generally not measurable directly, may be forth- 
coming from observations at widely separated frequencies, 
as emission at different frequencies is believed to origi- 
nate at different heights from the star leading to changes 
in beam size. For this, the dependence of the radiation 
frequency on the height, the so called radius-to-frequency 
mapping, should be known a priori. Alternatively, it is pos- 
sible to use the data on an ensemble of pulsars sampling a 
range of impact parameters. However, it is important that 
all the pulsars in the sample form a homogeneous set in 
terms of the profile types etc. Several attempts to model 
the pulsar beam have used the latter approach. Based on 
their study, Narayan and Vivekanand (1983) concluded 
that the beam is elongated in the latitude. Lyne & Manch- 
ester (1988), on the other hand, have argued that the beam 
is essentially circular (see also Gil & Han 1996, Arendt & 
Eilek 1998). Based on the dipole geometry of the cone of 
open field-lines. Biggs (1990) found that the beam shape 
is a function of the angle (a) between the rotation and the 
magnetic axes. The reasons that all these analyses predict 
different results could be manifold. For example, Narayan 
& Vivekanand used a data set consisting of only 16 pulsars 
and assessed the beam axial ratio on the basis of the to- 
tal change in the position angle of the linear polarization 
across the pulse profile. Apart from poor statistics, their 
analysis suffered from the large uncertainties in the polar- 
ization measurements available then. Lyne & Manchester 
(1988) used a much larger data set in comparison and 
examined the distribution of normalized impact parame- 
ter Pn = Pgo/pgo, where /Jgo & pgo are the impact angle 
and the beam radius computed for a = 90°. Based on 
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their observation that the distribution of /3„ is 'essentially 
uniform', they concluded that the beams are circular in 
shape. The apparent deficit at large /?„ is attributed to a 
luminosity bias. It is worth noting that the deficit is seen 
despite the fact that /3„ overestimates the true P/p (given 
that they disregarded the sign of /?), this is particularly 
so at large /3 values. 

Biggs (1990) used the same data set as well as the 
Pn distribution as used by Lyne and Manchester (1988), 
but drew attention to a 'peak' in the distribution at low 
Pn- The shapes of the polar cap defined by the region of 
open field lines, as derived by Biggs, show that the beam is 
circular for an aligned rotator, but undergoes compression 
along the latitudinal direction with increasing inclination 
a. 

In this paper, we address this question within the ba- 
sic framework advanced by Rankin (1993a) which, at the 
least, is qualitatively different from that of Lyne & Manch- 
ester (1988). The classification scheme (Rankin, 1983a), 
based on the phenomenology of pulse profiles, polarization 
and other fluctuation properties etc., provides a sound ba- 
sis for explicit distinction between the core and the conal 
components, with each of them following a predictable 
geometry (see also Oster & Sieber 1976; Gil & Krawczyk 
1996 for conal beams). Lyne & Manchester (1988), on the 
other hand, prefer to interpret the observed variety in 
pulse shape and other properties as a result of patchy il- 
lumination, rather than any particular pattern within the 
radiation cone. The observed differences in the properties 
of pulse components are then to be understood as grad- 
ual changes as a function of the distance from the center 
of the basic emission cone. Their analysis thus naturally 
disregards the possible existance of conal features. 

Assuming the possibly confined 'conal-componcnt' ge- 
ometry and by accounting for all the relevant geometrical 
effects, we re-examine the shape of pulsar beams and their 
frequency dependence. Recently published multifrequency 
polarization data, at six frequencies in the range between 
234-1642 MHz (Gould & Lyne, 1998), has made this in- 
vestigation possible. 

2. Data set 

For the present investigation requiring reliable estimates 
of a & /3, we use the data set comprised of only those 
pulsars whose pulse profiles are identified as 'triple' (T) 
or 'multiple' (M), as classified by Rankin (1993a, 1993b). 
The reason for the choice is that the T and M pulsars 
show a core component in addition to the conal compo- 
nents, so that a reliable estimation of the angle (a) be- 
tween the rotation axis and the magnetic axis is possible, 
using Rankin's (1990) method. In this method, the ratio of 
the observed core- width to the limiting width (2.45°P~°-^) 
is interpreted as the geometric factor 1/ sin(Q!), providing 
by far the most reliable estimates of a. For the conal dou- 
bles and conal singles, devoid of any core component, the 
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Fig. 1. Schematic representation showing the geometry of 
the pulsar emission region. 

estimates of a are less reliable. The core singles are natu- 
rally excluded from this analysis of the conal emission ge- 
ometry. For each pulsar in our selected sample, we define 
the conal width as the separation between the peaks of the 
outermost conal components. It is important to note that 
the nominally 'central' core component, which is argued 
to originate closer to the stellar surface, may not neces- 
sarily be along the cone axis. Such a possibility is clearly 
reflected in many pulse profiles where the core compo- 
nent is displaced from the 'center' definable from the conal 
components. Hence, the location of the core component 
is disregarded in our estimation of the conal separation. 
Columns 1 and 2 of table 1 list the name and profile type 
of these pulsars. Columns 3 to 8 list the calculated widths 
of the pulsars at frequencies 234, 408, 610, 925, 1400, and 
1642 MHz respectively. Column 9 gives the pulsar period 
in seconds. Columns 10 and 11 list the a and /? values of 
the pulsars taken from Rankin (1993b). 

Rankin (1990) has estimated the inclination angle a 
using the relation, sin(Q:) = 2.45°P~'^'^/Wcore) where 
Wcoro is the half-power width of the core component (at a 
reference frequency 1 GHz) and the period P is in seconds. 
The impact angle /3 has been estimated based on the ro- 
tating vector model of Radhakrishnan & Cooke (1969), 
using the relation sin(/3) = (dx/(i0)niax/ sin(a), where 
(rfx/c^<A)max is the maximum rate of change of the polar- 
ization angle x with respect to the longitude (/). 

In the following analysis, we treat the different fre- 
quency measurements on a given pulsar as 'independent' 
inputs much the same way as the data on different pul- 
sars, since the pulsar beam size is expected to evolve with 
frequency. Thus, at different frequencies one obtains inde- 
pendent cuts (at different /3/p) across the beam, though 
/3 remains constant for a given pulsar. This increases the 
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Fig. 2. Schematic representation of an elliptic pulsar 
beam of axial ratio R with the longitudinal and the lat- 
itudinal axis as h and Rb respectively. 5r is the width of 
the emission beam cone. See text for discussion on the 
connection between the 'gap-angle' 6g and 5r/r. 

number of independent constraints by a usefully large fac- 
tor. In fact, we would like to contrast this approach with 
the one where, for each pulsar, one obtains a best fit fre- 
quency dependence of the observed widths and then uses 
the data to obtain the width at a chosen reference fre- 
quency. The latter approach fails to take into account the 
dependence of the observed widths on (3/ p that is inherent 
for any non-rectangular shape of the beam. 

3. A direct test for the shape of beams 

The Fig is a schematic diagram illustrating the geome- 
try of pulsar emission cone. The emission cone, with half- 
opening angle sweeps across the observers line-of-sight 
with an impact parameter (distance of closest approach to 
the magnetic axis) (3. The spherical triangle PQS (refer to 
Fig. |l|) relates the angles a, (3 and the profile half-width 
(j^i, to the beam radius pv by the following relation (Gil, 
Gronkowski & Rudnicki 1984), 

sin^{pu/2) = sin2(0^/2) sin(a) sin(a + f3) + &\n^{fi/2) (1) 

The subscript v in p^, and (p^ denotes that these quantities 
depend on frequency v. This equation assumes that the 
cone is circular, in which case Pi, becomes independent of 
p. 

In reality, the beam may not be circular, but rather el- 
liptical with, say, R the axial ratio and b the longitudinal 
semi-axis of the ellipse as shown in Fig. ||. It is easy to see 
that the length of the radius vector r depends on the angle 
9 (with the longitudinal axis) when R is not equal to 1. 
The variation of r as a function of 6 for three different R 
values (namely 1, 1.5 and 0.5) are shown as examples in 
Fig. H. The p^, determined assuming that the cone shape 
is circular (as in Rankin 1993b) is indeed a measure of 
the radius vector r, once the period and frequency depen- 
dences are corrected for. Such data on (r,6') spanning a 
wide enough range in can therefore be examined to seek 
a consistent value of the axial-ratio R. However, if R is 
a function of a, as suggested by Biggs (1990), then the 




Fig. 3. The above curves illustrate the normalized vari- 
ation of r with 9 (refer to figure ^ with three different 
values of R. 



{t,9) samples would show a spread bounded by the curves 
corresponding to the maximum and minimum values of R. 

Such an examination of the present data suggests a 
spread below the line for i? = 1, indicating that the beam 
deviates from circularity and that the spread could be due 
to the a dependence of R. However, this deviation from 
circularity is not very significant. We discuss this in detail 
later in section 5. 

We have also examined the Pi, values obtained by 
Rankin (1993b) through such a test. However, no signif- 
icant deviation from circular beams was evident. We be- 
came aware of a similar study by C.-I. Bjornsson (1998), 
also with a similar conclusion. We note that the only dif- 
ference between our estimates of pu and those of Rankin 
is in the definition of the conal widths. Rankin defines 
the width as the distance between the outer half-power 
points (rather than the peaks) of the two conal outrid- 
ers, and the widths were then 'interpolated' to a refer- 
ence frequency of 1 GHz. Such estimates are prone to 
errors due to mode changes, differing component shapes 
etc., and to the effects of dispersion & scattering (some 
of which she attempted to accommodate). We measure 
the widths as the peak-to-peak separations of the outer 
conal components, which are less sensitive to the sources 
of error mentioned above. We have also confirmed (in the 
PSRs 0301-1-19, 0525-f21, 075H-32, 1133+16, 1737-^13, 
2122+13 and 2210+29 using the data from Blaskiewicz 
et al. 1991) that the 'peaks' of the conal components are 
symmetrically placed with respect to the "zero-longitude" 
(associated with the maximum rate of change of the po- 
sition angle), which is not always true for the outer half- 
power points. 
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Table 1. The table lists the pulsar name and the widths measured at 6 different frequencies from the observations 
of Gould & Lyne (1998). In several cases the widths could not be estimated due either to poor quality profiles or 
to absence of data. The a, p values are taken from Rankin (1990, 1993b). LM indicates that the (3 value (for PSR 
0656+14 and 1914+09) is taken from Lyne & Manchester (1988). 
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4. The model of the pulsar beam 

We model the pulsar beam shape as elliptical in general 
and express it analytically as, 

sin^((/>^/2)sin(a) sin(a + /3) sin^(/3/2) ^ ^ 

sin2(p,/2) sin2(i?p,/2) 



While a, (3 and (f)^ can be estimated, directly or indirectly, 
from observations, R and p,^ are the two parameters which 
in turn define the beam shape and size — and the available 
data set of T and M-profiles is expected to sample most of 
the I (3/ pv I range (0-1) with reasonable uniformity. The 
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imphcit assumption in this statistical approach is that a 
common description for i? & is valid for all pulsars. The 
common description should, however, account for relevant 
dependences on quantities, such as frequency, period, a, 
etc. properly. 

4-.1. Frequency dependence oj 

The radio emission at different frequencies is expected to 
originate at different altitudes above the stellar surface, 
with the higher frequency radiation associated with re- 
gions of lower altitude. This phenomenon known as radius- 
to-frequency mapping, finds overwhelming support from 
observations. Thorsett (1991) has suggested an empirical 
relation for the observed pulse width as a function of fre- 
quency, which seems to provide adequate description of 
the observed behaviour. We adopt a similar relation for 
the frequency evolution of the beam radius pi, as follows 

P, = p{l + Ku-^), (3) 

where p is the value of pi, at infinite frequency, C the spec- 
tral index, and K a constant. Note that both ( k K are 
expected to have positive values, so that the minimum 
value of pi, is p, which should correspond to the angular 
size of the polar cap. 

4-. 2. Period dependence on p^, 

Rankin (1993a) has demonstrated (see also Gil, Kijak & 
Seiradakis 1993; Kramer et al. 1994) that the beam radius 
p varies as p^^.s (.^^j^gj-g p jg i^i^g pgriod of the pulsar), a 
result which is in excellent agreement with that expected 
from a dipole geometry (Gil 1981). Eq ^ thus takes the 
form 

p,^Po{l + K,y-^)p-°-^ (4) 
where po is the minimum beam radius for P — 1 sec. 

4-.3. Functional dependence of R on a 

Biggs (1990) has suggested that R should be a function of 
a, such that the beam shape is circular for a — and is 
increasingly compressed in the latitudinal direction as a 
increases to 90°. We therefore model the functional depen- 
dence of i? on a as i? = i?oT, where Ro is the axial ratio 
of the beam at a = 0, and r is a function of a. According 
to Biggs (1990), i?o = 1 and r is given by 

T{a) = 1- Kix 10"^a - K2 x IQ-'^'a'^, (5) 

where Ki, K2 are constants and a is in degrees. Biggs 
finds that Ki and K2 are 3.3 and 4.4, respectively. We, 
however, treat Ki^2 as free parameters in our model. 



4-.4- The number of hollow cones 

Based on the study of conal components, Rankin (1993a) 
has argued for two nested hollow cones of emission- 
namely, the outer and the inner cone. Assuming the beams 
to be circular in shape, opening half angles of the two 
cones at 1 GHz were found to be 4.3° and 5.7°, respec- 
tively. 

During our preliminary examination of the present 
sample, we noticed a need to allow for three cones of emis- 
sion. To incorporate this feature in our model, we intro- 
duce two ratios, rl < 1 and r2 > 1, to define the size 
scaling of the inner-most and the outer-most cone, respec- 
tively, with reference to a 'middle' cone, for which the 
detailed shape is defined. 

Using the model here defined, we need to solve for Ro, 
C, K, Po, Ki, K2, rl and r2 in this three-conal-ring model. 
The parameter set thus represents an 'average' description 
of the beam. 

5. Results and Discussion 

An optimized grid search was performed for suitable 
ranges of the parameter values and in fine enough steps. 
For C,, the search range allowed for both -|-ve and -ve val- 
ues. By definition, ri < 1 and r2 > 1. The best fit was 
obtained by minimizing the standard deviation do defined 

by 

l Y.um ^. 180° 

where Di is the deviation of the data point from the 
nearest conal ring in the model and Ndof denotes the de- 
grees of freedom. The factor ISO/tt gives CTo in units of de- 
grees under the small-angle approximation. Table 2 lists 
the parameter values which correspond to the best fit for 
the entire sample set. With these values, the eq. || can now 
be rewritten as 

p. = 4.8°(l-H66i.MHz)^"°', (7) 

where pi, is in degrees. This average description for the 
'middle' cone applies also to the other two cones when 
Pi, is scaled by the ratio rl = 0.8 or r2 = 1.3 (for the 
inner and the outermost, respectively). Fig ^ shows the 
data (plotted to a common scale) for one quadrant of the 
beam and the three solid curves corresponding to the best 
fit cones. The points in the figure, though corresponding 
to different pulsars and frequencies, are translated to a 
common reference scale appropriate for P = 1 sec, a = 
and i> = 00. 

We have assumed the period dependence of pi, as 
P^"-^, whereas Lyne and Manchester (1988) found a de- 
pendence of P~ 3 . We have examined the latter possibility 
and found that the difference in the standard deviation is 
at the level of 2.5-3 a and we cannot rule out the P^3 law 
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Table 2. The best-fit model parameters for the shape of conal beams. The error bars correspond to a la uncertainty. 
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ONE QUADRANT OF THE BEAM 




Fig. 4. Distribution of the (x,y) locations of the conal 
components on a common scale. The three solid lines indi- 
cate the three emission cones in the quadrant shown. The 
circles with crosses refers to pulsars with a values less than 
45° and the filled circles with a greater than 45°. 

Distribution of effective fi/p 



Fig. 5. Histogram of the distribution of effective ^. 



with confidence. We have also checked for the dependence 
of i? on a by using 3 sub-sets, each of range 30° in a. 
The best fit values for R in the different a segments are 
l±g:|, 0.8±g:| & 0.5±g:| for a ranges 0° - 30°, 30° - 60° 
& 60° — 90°, respectively. This dependence of Ro on a, 
even if it were significant, is quite consistent with our val- 



ues of Kl, K2 (Table 2) as well as with the results of 
Biggs (1990). However, given the uncertainties in the R 
estimates for the three ranges, it is not possible presently 
to rule out a dependence of R on a. Indeed, this part of the 
goodness-of-fit is negligible, do (the standard deviation) is 
0.18° when Kl and K2 ^ and 0.2° when K1,K2 = 0. 
Earlier Narayan & Vivekanand (1983) had argued that R 
is a function of the pulsar period. To assess this claim, 
our sample was divided into three period ranges and the 
corresponding R estimates compared. However, no period 
dependence was evident and it was possible to rule out 
such a dependence with high confidence. 

The number and thickness of conal rings: As 
already noted and can be seen in Figure ^, we do see 
evidence for a possible cone outside the two cones dis- 
cussed by Rankin (1993a). Also, presence of a 'further 
inner' cone has been suggested by Rankin & Rathnasree 
(1997) in the case of PSR 1929+10. The pulsars sugges- 
tive of this outer cone (refer Figure]^) are PSRs 0656+14, 
1821+05, 1944+17 and 1952+29 (at frequencies 234 MHz 
and higher). We have examined the possibility that these 
cases really belong to the central-cone, but are well out- 
side of it due to an error in the assumed values of a. We 
rule out the possibility as the implied error in a turns 
out to be too high to be likely. It is important to point 
out that a noisy sample like the present one would ap- 
pear increasingly consistent, judging by the best-fit crite- 
rion, with models that include more cones. The question, 
therefore is whether we can constrain the number of cones 
by some independent method. In this context, we wish to 
discuss the noticeable deficit of points at high P/po- Since 
the deficit reflects the absence of conal singles and conal 
doubles in our data set, the size of the related 'gap' at 
large 9 values, can be used to estimate the possible thick- 
ness of the conal rings. The absence of points at 9> 60° 
(Figure ^) suggests that the conal rings are rather thin, 
since a radial thickness Sr comparable to the ring radius 
would imply a wider gap in 9. To quantify this, we write 
the following relation, 



Sr — 2r 



(l-sin0,) 



(1 



smf 



,)' 



(8) 



where 9n is the 9 at the start of the gap (as illustrated 
in Fig. |). With 9g - 60°, Sr/r would be about 20%. 
The presence of more than one distinguishable peak in the 
distribution of beam radii (shown in the bottom panel of 
Fig. ^ clearly indicates that the conal separation is larger 
than the cone width. This combined with our cone-width 
estimate suggests the number of cones is 3 (for the present 
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range of radii), providing an independent support for our 
model. This picture is consistent with the estimates by Gil 
& Krawczyk (1997) and Gil & Cheng (1999). 

Component separation vs. frequency: It is inter- 
esting to note that for certain pulsars the cone associated 
with the emission seems to change with frequency. For 
example, the conal emission in PSR 1920+21 appears to 
have 'switched' at 610 MHz to the innermost cone while 
being associated with the central cone at other frequen- 
cies. Rankin (1983b), in a comprehensive study of the de- 
pendence of component separation with frequency, invokes 
deep 'absorption' features to explain the apparent anoma- 
lous reduction in the component separation in certain fre- 
quency ranges. We suggest that such anomalous reduction 
in the separations could be due to switching of the emis- 
sion to an inner cone at some frequencies. Observations 
at finely spaced frequencies in the relevant ranges would 
be helpful to study this effect in detail. The other pul- 
sars which show similar trends are PSRs 1804-08, 2003- 
08, 1944-hl7 and 1831-04. It should be noted that such 
switching is possibly reflected, also, in mode changes. 

The deficit at low (3/ po- The absence of points near 
/3 = is clearly noticeable in Fig. ^. Such a 'gap' is also 
apparent in the distribution of /3/po plotted in Fig. ||. 
The gap was already noted by Lyne & Manchester (1988). 
They argued that it arises because the rapid position-angle 
swings (expected at small /3's) are difficult to resolve due 
to intrinsic or instrumental smearing, leading to underes- 
timation of the sweep-rates. With the improved quality 
of data now available, the intrinsic smearing is likely to 
be the dominant cause for this circumstance. There are a 
number of clear instances among the general population of 
pulsars where the polarization angle traverse near the cen- 
tral core component is distorted. PSR 1237-1-25 provides 
an extreme examples of such distortion, and Ramachan- 
dran & Deshpande (1997) report promising initial efforts 
to model its polarization-angle track as distorted by by a 
I0W-7 core-beam. Another possibility for the low-/3/po gap 
is that it could simply be a selection effect caused by less 
intense emission in the cone center than at intermediate 
traverses. If so, the low frequency turn-overs in the energy 
spectra of pulsars may at least be partly due to this, since 
at lower radio frequencies the (3 / po is relatively smaller. 

The sources of uncertainties in the present anal- 
ysis: The standard deviation cto corresponding to the 
best-fit model amounts to about 15% of the conal radius. 
This fractional deviation (comparable to the thickness of 
the cone) is too large to allow any more detailed descrip- 
tion of the beam shape (such as dependence on a, for ex- 
ample). We find it useful to assess and quantify the sources 
of error, partly to help possible refinement for future in- 
vestigations. The three data inputs to our analysis are 
a, (3 and (jii,, while the basic observables are the maximum 
polarization-angle sweep rate and core width, apart from 
the conal separation measured. It is easy to see that the 
errors in the core-widths will affect directly both a and 



P estimates. Over the range of 9s spanned by the present 
data set the errors in a are likely to dominate, since the x 
& y (in figure 1^) are almost hnearly proportional to sin(a). 
Hence, the fractional deviation may be nearly equal to (or 
define the upper limit of) the fractional error in sin(a) 
and therefore in the core-width estimates. 

Rankin (1990, 1993b) notes that in several cases the 
apparent core-widths might suffer from 'absorption' and 
the widths might be underestimated if the effect is not 
properly accounted for. Also, in some cases, the widths 
were extrapolated to a reference frequency of 1 GHz using 
a dependence. There have been several suggestions 

regarding the 'appropriate' frequency dependence which 
would give significantly different answers when used for 
width extrapolation. For example, if our best-fit depen- 
dence for conal width is used for the core-width extrap- 
olation, the values would differ from Rankin's estimates 
(through extrapolation) by as much as 15%, enough to 
explain the present deviation in some cases. Another pos- 
sible source of error is the uncertainty in the sign of (3 
(important only for the sin(Q! + f3) term in equation and 
hence for small a). As Rankin points out, it is difficult 
to determine the sign unambiguously in most cases and 
hence the information is only available for a handful of 
pulsars. 

Evidence in favour of 'conal' emission: The signif- 
icant implication of the gap at 9> 60° (referred to earlier) 
deserves further discussion. If the 'conal' components were 
results of a merely patchy (random) illumination across 
the beam area, (as argued by Lyne & Manchester, 1988), 
then such a gap should not exist. If a single thick hollow 
cone were to be responsible for the conal components, a 
gap (corresponding to the conal-single types) would still 
be apparent but then it should be above a cut-off y value 
(refer figure 4) and not in a angular sector like that ob- 
served. On the other hand, if indeed the conal emission 
exists in the form of nested cones (as distinct from the 
core emission), then the shape of the gap is a natural con- 
sequence of our not including conal-single profiles in this 
analysis. This gap, therefore, should be treated as an im- 
portant evidence for a pulsar beam form comprised, in 
general, of nested cones of emission. 

6. Summary 

Using the multifrequency pulse profiles of a large number 
of conal-triple and multiple pulsars we modelled the pulsar 
beam shape in an improved way. Our analysis benefits 
from the different frequency measurements being treated 
as independent samples, thus increasing the number of 
independent constrains. The main results are summarized 
below. 

1) Our profile sample is consistent with a beam shape 
that is a function of a, circular at a = and increasingly 
compressed in the latitudinal direction as a increases, as 
suggested by Biggs (1990). However, the data is equally 
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consistent with the possibility that the beam is circular 
for all values of a. 

2) We identify three nested cones of emission based 
on a normalized distribution of outer components. The 
observed gap {6> 60°) in the distribution independently 
suggests three cones in the form of annular rings whose 
widths are typically about 20% of the cone radii. We con- 
sider this circumstance as an important evidence for the 
nested-cone structure. 

Any further significant progress in such modelling 
would necessarily need refined estimates of the observ- 
ables, particularly the core-widths. 
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